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Abstract 

In this paper we compute the most general nondiagonal reflection matrices of 
the RSOS/SOS models and the hard hexagon model using the boundary Yang- 
Baxter equations. We find new one-parameter family of the reflection matrices for 
the RSOS model in addition to the previous result obtained in ||. We also find 
three classes of the reflection matrices for the SOS model, which has one or two free 
parameters. For the hard-hexagon model which can be mapped to RSOS(5) model 
by folding four RSOS heights into two, the solutions can be obtained similarly with 
a main difference in the boundary unitarity conditions. Due to this, the reflection 
matrices can have two free parameters. We show that these extra terms can be 
identified with the 'decorated' solutions. We also generalize the hard hexagon 
model by 'folding' the RSOS heights of the general RSOS(p) model and show that 
they satisfy the integrability conditions such as the Yang-Baxter and boundary 
Yang-Baxter equations. These models can be solved using the results for the RSOS 
models. 
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1 Introduction 



In the study of the two-dimensional integrable models of quantum field theories and 
statistical models, the Yang-Baxter equation (YBE) plays essential roles in estabilish- 
ing the integrability and solving the models. Recently there has been a lot of efforts in 
introducing the boundaries into the integrable systems for possible applications to con- 
densed matter physics and statistical systems with non-periodic boundary conditions. 
The boundaries entail new physical quantity called reflection matrices which depend 
on the boundary properties. 

The boundary Yang-Baxter equation (BYBE) (also known as the reflection equa- 
tion) |Q] is the necessary condition for the integrable statistical models |2|, ^] and quan- 
tum field theories Q with boundary. The equation takes the form 

Ri(u)Su(u + u)R2(u)Si 2 (u -u) = Sn(u -u)R 2 (u)Sn(u +u)Ri(u) (1.1) 

where Ri(2) 1S the boundary reflection matrix in the auxiliary space 1(2) and £12 is the 
solution to the YBE. 

To date, several solutions of the BYBE have appeared in the literature. Compared 
with the vertex-type models, however, far less is known for the solution in the face- 
type model such as the solid-on-solid (SOS) or restricted-solid-on-solid (RSOS) model. 
Among these, most of the known solutions are 'diagonal' in the sense that the reflection 
matrices are diagonal|5|, ||. Using these, one can find the partition functions in the 
infinite lattice limit 0. 

Nondiagonal reflection matrices are more interesting because they can be applied to 
physical phenomena. Furthermore, one can find explicit vertex-face correspondence if 
one can classify nondiagonal solutions completely. In addition to those found in Q , we 
will show that different classes of nondiagonal reflection matrices are possible. These 
new solutions include one or two free parameters which are related to those in the 
boundary potential. The explicit formulae between two sets of parameters are still not 
clear. We hope our complete nondiagonal solutions for the SOS model may be related 
explicitly with those of the boundary sine-Gordon model. 

^From a physical point of view, many interesting models are face- type; the RSOS/SOS 
models, the hard-hexagon model (HHM) etc. These models play very important roles in 
the statistical mechanics system and in the quantum field theories such as the perturbed 
conformal field theories. 

In this paper we derive complete nondiagonal reflection matrices for the RSOS/SOS 
and the HHM in a unified way. We will express BYBE as a linear equation which should 
satisfy extra nontrivial conditions. Due to the linearity, the general solutions are linear 
combinations of each solution with arbitrary coefficients. Some of these coefficients are 
fixed by the boundary crossing and unitarity conditions The reflection matrices 
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of the HHM can be constructed from those of the RSOS model since the HHM can 
be mapped to the RSOS (5). We consider a similar mapping for the generic RSOS(p) 
model and their boundary reflection matrices. 



2 RSOS(p) model 

In this section we solve the BYBE for the RSOS(p) ; p = 3, 4 . . . scattering theory. The 
RSOS(p) scattering theory is based on a (p — 1) - fold degenerate vacuum structure, 
vacua can be associated with nodes of the A p -i Dynkin diagram. The quasi particles in 
the scattering theory are kinks that interpolate neighboring vacua, they can be denoted 
by non-commutative symbols K ab (u) where \a — b\ = 1 with a,b = 1, . . . ,p — 1 and u 
is related to the the kink rapidity 9 by u = —i6/p, so that the physical strip is given 
by < Re u < tt /p. In the rest of the paper, we will refer to a, b as heights or spins. 
Formally, scattering between two kinks can be represented by the following equation 

K da (u)K ab (u) = ]T Sf c {u - u)K dc {u)K cb (u) , (2.1) 

c 

where the S'-matrix is given by 



S£{u) = U{u) ( |m ) U/27 W£(v) . (2.2) 



The Boltzmann weight, 

Wt (u) = smu5 bd (j*) ' + sin( 7 - u)5 ac , (2.3) 



satisfies the YBE in the RSOS representation. 
Here [a] denotes the usual g-number, 

sin(<27) 7r 

[a] = — 7 = - , 

sin 7 p 

and the overall factor U{u) is a product of Gamma functions, 

mu) = Ir^rfi-^rfi-l + ^n^"'^ 7 ""' 



TT \ TT I \ TT 



TT TT J ~_~ Ff(0)Ff(7) ' 



F,( u ) - ^ ^ ^ - ^ (2 4) 

il J r f (21+1)7 uWj ■ (2i-l)7 u\ ■ 

\7T 7r /\ 7r / 
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This factor, satisfying the relations 



U(u)U(— u) sin(7 — u) sin(7 + u) = 1 

W(7 - u) = U{u) , 

together with the overall (/-number factor ensures that the 5-matrix satisfies both cross- 
ing and unitarity constraints: 

S£( 7 -u) = S£{u) (2.5) 
Y,Sf c ,(u)sl\-u) = 5 ac . (2.6) 

c 

Let us consider now the above scattering theory in the presence of a boundary. The 
scattering between the kink and the boundary denoted formally by B a is described by 
the equation 

K ab (u)B a = £ R b ac (u)K bc (-u)B c . (2.7) 

c 

Note that in this representation, the boundary naturally carries an RSOS spin. 

The function R b ac is called the boundary reflection matrix and satisfies the BYBE, 
which in the RSOS representation takes the form 

£ R a bb , (u)S a / a , (u + u)B$ b „ («')S* «„ (u -u) = 

a ,b 

^ S£ (u - u)Ri, (u')S°% (u + u)Ri' b „ (u) . (2.8) 

a ,b 

In general, the function R bc (u) can be written as 

R a bc (u) = K{u) [j^) [h^cXUu) + 4c {S b , a+ iU a (u) + <5 6)0 _iD„(«)}] , (2.9) 

where 7?.(n) has to be determined from the boundary crossing and unitarity constraints, 
while X bc and U a , D a have to be determined from the BYBE. An overall g-number factor 
has also been multiplied to the above to cancel that from the bulk 5-matrix in order 
to simplify the BYBE. If X bc does not vanish, the boundary i?-matrix describes non- 
diagonal scattering process, otherwise the scattering is called diagonal. Note that due 
to the restriction that vacuum assumes value 1, ... ,p — 1, X bc , -Xf" 1 , D±, U p -i are not 
defined. The case p = 3 has only diagonal scattering, so X bc does not exist. 
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We will concentrate in what follows to the scattering where the off-diagonal com- 
ponent Xg c is non- vanishing. To start, the case b / c ^ b in eqn. fl2.8Q gives 

K-i,a+i(u')X a a +l a+3 (u) = X-_ 1:a+1 (u)X a a +l a+3 (u) ; 2 < a < p - 4 . (2.10) 

This equation implies that X£±i can be written as 

X a a±1 , aTl (u) = h ± (u)X a ± , 

where h±(u) depends only on u and X^_ only on a. 

On the other hand, the case c = b = b" , a = a" gives 

X a a -i,a + i{u )K+i,a-i{u) = ^ W W^i,„-i(«') ; 2 < a < p - 2 , (2.11) 
which implies that 

h + (u )/i_(u) = h + {u)h-{u ) , 
from which we conclude that 

h + (u) = (const. )h^(u) . 

Absorbing the constant in the above equation into X a or JT", we can make h + equal 
to /i_ so that we can absorb the h±{u) into the overall 7Z(u) factor and treat X£ c as u 
independent from now on. 

With this simplification, eqn.(|2,8j) can be broken down into the following indepen- 
dent equations in addition to the above two equations: 

Kd^u (i + /-pfeLj) + d: m . (i + /+ £±! 

+ X a+l,a+3 X a+i,a+lf- = U a D a+2 f+ + f- r + ^ 



+ C7- a Z7 a /_ ( 1 + /+ r^T ) + K-i,a + iX a a+ i,a-if- (2-12) 



for 1 < a < p — 3, 



w-(i+/ + I ^)+f:- 1 / + fi+/- 1 ^ 



aj 

E/a-i/+ - CWi/- (2-13) 



r„/-(i-i-/ + r^TT]+A (+ >/+('+./ + 



[a + l]7 a+ ^ + V [a + 1], 

D a+2 f + - D a f^ (2.14) 
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for 2 < a < p — 3, and 



[o-l] 2 " "V [a-l]A "[o-l] 



[o+l] 2 " "\ '~[a+l\) V [a + l] 

for 2 < a < p — 2. In the above equations, we used a compact notation where U a = 
U a (u), U a = U a (u) (similarly for D a ) and 

f± = sin(u ±u)/ sin(7 — u =F u ) • 

In addition, it should also be mentioned that the last term in the rhs (lhs) of 
eqn.(|2.12 ) is present only when a ^ l(p — 3) and the first terms of eqns,( ^T5| ) and 
( (2.16Q are allowed only for a ^ 2 and a ^ p — 2, respectively. 

Before solving the equations directly, it is helpful to investgate the structure of 
BYBE. First of all, the BYBE is covariant under the transformation 

a -> p - a , U a (u) ->• ±D p _ a (u) , (2.17) 

and have the following symmetry 

U a {u) = -D a {-u) (2.18) 

for 2 < a < p — 2. Futhermore, since the eqns.( ^13 )-( 2.1^ ) are linear, their general 
solutions are linear combinations of 'fundamental' ones. Another important fact is that 
the amplitudes with a even and those with a odd are completely decoupled in BYBE 
and give different solutions in general. However if p is odd, two sets are related by 
( p. 17 ) and have the same solutions. 

By solving the linear equations ( 2.13[) -( p.l6 ), we find the most general nondiagonal 
solution for p > 5: 

U a (u) = A sin(2n + 07) + + 

sm 2u sin 07 

| ep-l£ f sin(2n + 07) _ 

sina7 \ sin2u 
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D a (u) 



A sin(2u — 07) + 



B _ e p e a -iC 
sin 2u sin a 7 



(2.19) 



e p -iD ( sin(2u — 07) 



sin 07 



sin 2u 



-1) 



where e a is 0(1) if a is odd(even). and A, B, C and D are free parameters. 

Having found U a , D a , the function X£ c can be easily obtained from eqn. ( EU2D , after 
taking u to be — it since X^, does not depend on the rapidity. This gives 

Ktla+sKtla+l ~ K-l,a + lK+l,a-l = U a (-u)U a {u) - D a+2 {-u)D a+2 {u) 

for 1 < a < p — 3. Substituting U a , D a into the rhs and iterating the equations, we get 

K-i,a+iK+l,a-i = e P e a Xf 3 Xi 1 + A 2 {sin 2 (l + e p e a ) 7 - sin 2 07} 

—2AB {cos(l + e P e a )7 — cos 07} 



~2 

+£p£a-l^ 



1 



1 



sin 2 7 sin 2 07 



(2.20) 



+2e p ^D 2 



sin 2 7 sin 2 07 



+2ep_i£> 



2 J cos 7 



sin 2 7 



+ (-1) 



a cosa7 
sin 2 07 



Since this equation fixes only the product, X^_ x a+1 and X^ +1 a _i are determined upto 
a gauge factor. Inserting a = 1 one gets X^X^q = as expected. For p even, Xj^Xfi 
is not determined yet. 

Consider now boundary unitarity and crossing symmetry conditions for the reflec- 
tion matrix R^ c (u). Due to these conditions, the overall factor R(u) should satisfy 



= s bd 

c 

£S£(2u)/&( 7 /2 + u) = i? 6 a c ( 7 /2-«). 
d 

In terms of (|2.9| ) and (2.18), the unitarity condition becomes 

K{u)n{-u) [x a a+ha „ x X^ a+1 - U a (u)D a (u)] = 1 ; 2 < a < p - 2 
K^Ki-ujUi^Uxi-u) = 1 
^(u)^(-u)L> p _i(n)Z)p_i(-n) = 1 . 



(2.21) 
(2.22) 
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By inserting eqns. plS| ) and p3c| ) to the above, we find the following nontrivial con- 
straints for the free parameters: 



AB = AC = AD = 

C 2 
sin 2 7 



ii) Xf^X^! = A 2 (sin 2 7 — sin 2 2j) H — — — for even p . 



Note that all the cases satisfying the constraints, have at most one free parameter since 
we can absorb overall constant into lZ(u). We list each class of the reflection matrices 
as follows: 

Class (I) for general p (p ^ 3, 4): B = C = D = 0, A=l 

U a (u) = sin(2u + cr/) 

D a (u) = sin(2n-a7) (2.23) 
K-i,a+i X a+i,a-i = sin 2 7 - sin 2 a 7 . 

These weights have U p - a (u) = —D a (u) symmetry and no free parameter. This 
solution is the one obtained in ||. The unitarity condition gives 

n{u)TZ{-u){- sin 2 2u + sin 2 7) = 1 . (2.24) 

The crossing symmetry condition becomes 

U(2u)K(-f/2 + u) sin( 7 - 2u) = K(j/2 - u) . (2.25) 

The factor lZ(u) can be determined from eqns. ( ^24| ), ( |2.25| ) up to the usual CDD 
ambiguity by separating TZ{u) = 1Zq{u)1Zi(u) where TZq satisfies 

K (u)K (-u) = 1 
U(2u)K {7/2 + u)sm(j -2u) = Kq^/2 - u) , (2.26) 

whose minimal solution reads 

F (u) 



Ko(u) 



F (-u) 
While 1Z\ satisfies 



^i(u)7ei(-n)(-sin 2 2n + sin 2 7) = 1 

Ki(u) =Ki(j-u) (2.27) 
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with minimal solution 

^1 («) = ^(7/2, u)<j(7r/2 - 7/2, «) . 
Here cr(x, u) is a well-known building block satisfying the relations 

a(x,u) = a(x, r y — u) 
a(x,u)a(x, —u) = [cos(x + u) cos(x — u)]^ 1 , 



and is given by 
a(x, u) 



u(x, 1 - u) n(-x, 1 - u) n(x, -1 + u) n(-x, - j + U ) 
n 2 (x,i)n 2 (-x,i) 



Wx,u) = TT— 7^ ^ - x 



Class (II) for even p (p ^ 4): A = 0, C = 1 



sin 2u 


+ — — 


sin 07 




e a -i 


sin2n 


sin 07 


1 


e a -i 


sin 2 7 


sin 2 07 



D a (u) = — — -4^- (2.28) 



-^■0-1,0+1^-0+1,0-1 ■ 

This solution satisfies U p - a (u) = U p (u) (similarly for D a (u) and X?) and include 
one free parameter. To fix the overall factor TZ(u), 1Zo(u) is the same as Class (I) while 
the IZi (u) satisfies 

TZ^TZ^-u)!—, =1. (2.29) 

y sin 7 sin 2u J 

The minimal solution is 

. a{x,u)a(ir/2 - x,u) 

K^u) = sm 7 , 

cr(0, M)cr(7r/2, u) 

where 

sin2x = i?sin7 . 
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Class (III) for odd p (p + 1 3): A = 0, D = l 



u a = -^ + -^-( sm(2 \ +a7) -(-ir 

sin 2it sma7 L sin 2u 

D. = -+ --i_( si "' 2 " -^ _(_!).} (2.30) 
smzii sma7 L sin 2m J 

J 1 11 J cos 7 cos 07 

^a-l,a+l^a+l,a-l = ^ < — f + 2 S "T— § h ( — 1) "T— § 

L sin 7 sin aj J I sin 7 sin 07 

This solution has U p - a (u) = D p (v) symmetry and one free parameter. While Hq(u) 
does not change, the 1Z\{u) satisfies 



— ^ - — r-jr- = 1. 

\ sin tt sin 2tt sin 2tt / 



1 2B cos 2u 1 + 5 

3 

2 

The minimal solution is 

.7 cr(xi,n)<T(x2,n) 

/VI (it) =Sin— — r ; ; , 

where 

cos 2 x\ + cos 2 X2 = 1 + B sin 2 ^ 

1 7 
cos xi cos X2 = -(1 + 5) sin— . 

In the above analysis, we omit special cases of p = 3, 4 since p = 3 has only the 
diagonal reflection matrices and p = 4 has been extensively studied in ||. 



3 the SOS model 

We have considered in the beginning that the heights take values from 1 to p — 1, which 
is necessary for the bulk scattering weights to be finite as the parameter tt/^ = p is a 
positive integer. When tt/^ is not a rational number, there is no bounds on the heights 
and the corresponding representation is known as solid-on-solid (SOS). The removal 
of the heights' restriction certainly affects to set e p = 1 in the solutions ( |2.19| ) of the 
BYBE. Thus the solutions of BYBE in the SOS representation are 

B € C 

UJu) = Asin(2n + a7)H h a-1 

sin 2u sin 



9 



D f sin(2-u + 07) 



sin 07 [ sin 2k 

B e„_iC 



D a (u) = Asm(2u — aj) + — 



sin 2u sin 07 
I? [ sin(2ii — 07) 



sin 07 [ sin2u 

K-i,a+iK+i,a-i = e a X 2 3 X 2 1 +A 2 {sm 2 (l + e a h-sm 2 a 1 } (3.1) 

—2AB {cos(l + e a )7 — cos 07} 

+ £a-lC 2 (-T"2 ! o l 

L sin 7 sin 07 J 



-e a -iCT> 



+L> 2 



' 7 sin 07 . 
1 1 



cos 2 Z COS 2 ^? 



cos 2 (li^h COS 2 f 



redefining C as C — 2D. 

Inserting above solution to the unitarity condition restricts the coefficients in the 
same way as RSOS(p). We classify the solutions into three classes: 

Class (I) C = 0, A = 1 

j . , B D fsin(2u + a7) 

U a (u) = sm(2u + a 7 + — — + \ ^ V -1 

smiu sma7 L sm 2u 

D M = si n(2 u - a7 ) + -g r --^{ S "" 2 ";^ -l^ (3.2) 

sm 2u sma7 L sm 2u 

X a-i,a+i X a+i,a-i = sin 2 7 sin 2 C17 2i? (cos 7 cos 07) 

+W 1 1 ^ 



COS 2 3 COS 2 ^? 



The overall factor 7^o is the same as that of RSOS(p), but TZ\(u) now contains all 
the information of the boundary conditions and has to satisfy 

Hi (u)TZi(-u) - sin 2 2u - 2D cos 2u + sin 2 7 - 2B cos 7 + 



cos 2 ^ 



2BDcos2u B 2 + D 2 



: j , — , = 1 (3.3) 

$m z 2u siir2ti / 



K 1 (u)=1l 1 (i-u) 
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whose minimal solution is 



fti(u) 



2a(0,u)a(ir/2,u) 

where X1-X4 are related to B, D via 

4 

^cos2xj = —2D 



i=i 



cos2xjCos2xj 

i>j=l 



-2 + sin 2 7 - 25 cos 7 + 



cos 2 J 



cos 2xi cos 2xj cos 2xk = 2(1 + B)D 

i>jf>fc=l 

ry 

cos 2xi cos 2x2 cos 2x3 cos 2x4 = (cos 7 + B) 2 + D 2 tan 2 — 



Class (II) A = 0, C = 1 

A»(«) 



^0-1,0+1^0+1,0-1 



B gg-i -D f sin(2tt + 07) 
sin2u sin 07 sin 07 \ sin2u 

e a _i D f sin(2ii — 07) 



B 



sin2ii sina7 sin 07 t sin2n 
{D{1 - cos 7 ) - l} 2 {D(l - cosa 7 ) - e a _i} 2 



sin 2 7 



sin 07 



Now 

1li(u)1li(-u) 
whose minimal solution is 



{D(l - cos 7) - l} 2 ^ 2 2ff£> cos 2n B 2 + D 2 
sin 2 7 sin 2 2u sin 2 2m 



= 1 



sin 7 



a(xi,u)a(x 2 ,u) 
{D(l-cos 7 )-l} 2 + J D 2 sin 2 7 ff(0,«)<7(7r/2,«) ' 

where xi,X2 are related to B,D via 

(£ + £>) sin 7 

COS X\ COS X2 = 



COS 2 Xl + COS 2 X2 = 1 + 



2y/{D(l - cos 7) - l} 2 + D 2 sin 2 7 
i?D sin 7 



{D(l - cos 7) - l} 2 + Z? 2 sin 2 7 
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Class (III) A = C = 0, D = 1 



UJu) = » + ^{-202 + £2U 
smZu smo7 ^ smzu 

. . B If sin(2u — 07) , 

D a {u) = — : \ V- - 1 3.5 

smzu sma7 ^ smzu 



-^-a-l.a+l^a+l.a-l 



1 1 



Now 



Ki(u)lli(-u) 
with minimal solution 



cos 2 J cos 2 ^ 



1 2B cos 2m 1 + B 1 



cos 2 3 sin 2 2u sin 2 2u 



where 



7 o-(xi,u)cr(x 2 ,u) 

/Cl(u)=COS— — ; r ; ; r- 

u ; 2 ff(0,«)ff(7r/2,«) 



cos 2 X\ + cos 2 = 1 + B cos 2 ~ 

1 7 
cosxicosx2 = ~ (1 + -B) cos — . 



Note that the number of free parameters in the boundary reflection matrices of both 
vertex (the sine-Gordon model) and SOS Class (I), (II) representations are the same; 
two for non-diagonal and one for diagonal || ||. This strongly suggests that there can 
exist a well-defined transformation between the two models even with a boundary. 



4 Hard Hexagon Model 



The paticle spectrum of the HHM consists of a triplet of fundamental kink states 
Kqi,Kiq and KooP]- The bulk S'-matrix is given by [10, O] 



sfM = u(6) ( Ez 



\PdPb 



w£{u) 



(4.1) 
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with Boltzmann weights 

1/2 



*E^ + fc , (4,2) 



where 



p = 2 cos p, pi = 1 

7T 9t 
p = -,u=-0 

U(e)=Uo(6)Ui{u) 

sinh# — zsin| sinh # + i sin ^ 
sinh0 + isin| sinh (9 - i sin ^ 
sin(p — u) sin(2p + u) 

Ul(u) = —— r . tt r . 

sin(p + u) sm(2p — u) 

By mapping a = 2,3(1,4) of the RSOS(5) to a = 0(1) of the HHM, one can reproduce 
the bulk 5-matrix of the HHM from that of the RSOS(5). This means the RSOS(5) is 
homeomorphic to the HHM with differences in the overall factor U(u) and the relations 
between the spectral parameter u and the rapidity 9. 

These mean that two BYBE can be mapped to each other and the solutions of the 
HHM can be obtained by that of the RSOS. Writing the reflection amplitude of HHM 
as 



R a b M = n{9)( 



PbPc 



-e/2wi 



\PaPa, 

x [5 b ^ c (u) + 5 hc {6 bta+1 U(u) + 5 bta V(u) + <S 6)0 _i£>(u)}] , (4.3) 



the non-diagonal solution is 



, B D fsin(2u + 3p) 

U(u) = A S m(2u + 3p) + —— + ——{ \ W +l 

sinzti sm3p L smin 

r » , . B D fsin(2u — 3p) 

V(u) = A S m(2u-3p) + — — - \ 9 W +l 

sm 2u sm3p L sin 2u 

. „ . . B D fsin(2n + p) 1 

D(u) = Asm(2u + p) + —— + - — V . - — + l\ 4.4 

sin2w snip L sin 2m J 



X^X^ = A 2 (sin 2 n - sin 2 3p) - 2AB (cos p - cos 3p) 



-0 ^0 

— ^ ^ Sili , 

1 1 



+D 2 



sin 2 ^ sin 2 ^ 
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which can be easily read from eqns.( |2~T9| ),( |2.20 ). 



The unitarity and crossing symmetry conditions now reduce to 

K(9)K{-9)D(u)D(-u) = 1 (4.5) 

K(ni/2 -9) = U (26 f ni ^ + l U \ ll(iTi/2 + 9) . (4.6) 
sm(2^ — 2u) 

It is remarkable that for the HHM unitarity condition can not further reduce the 
arbitrary coefficients A, B and D. Among these the B and D terms are 'decorated' 
solutions which can be constructed from a fundamental solution R% c by 

R/d(«)V;cd = E S l> ~ ui)S a c e d {u + u x )Rt c {u). (4.7) 

a 

It is easy to check this satisfies the BYBE if S[^(u) is the solution of the bulk YBE 
with arbitrary u\. Using the trivial solution R^ c oc 5^ c , one can check that B and D 
terms can be obtained in this way. We will set therefore B = D = 0, A = 1 from now 
on. 



The overall factor 1Z(9) is can be determined from eqns.(pL5|),(4.6). Let 



such that 



U{9)=TIq{6)T1 1 {u) (4.8) 

K (9)K (-9) = 1, K {ni/2 -9)= U (29)n (Tri/2 + 9) 

Tlx (u)fti (-«)(- sin 2 2u + sin 2 /i) = 1 (4.9) 

„ , i \ sin(2u + 2u) „ , , 

Kl W2 -7T-U)= ■ / - - 7T + U , 

sm(2/i — 2u) 



then the minimal solutions are 

K (9) : 



sinh(| + f ) sinh(| - §) sinh(f + jg) sinh(| + g) sinh(| + lg) 



sinh(| - f ) sinh(§ + §) sinh(| - ^) sinh(| - §) sinh( 2 j( , , 
1 ; sin(/x + 2u) 

By generalizing the mapping of the RSOS(5) to the HHM, we can construct some 
generalized HHM whose particle spectrum consists of kinks K a b where \a — b\ = 1 with 
a, b = 0, . . . ,n — 1 and Kq$. Referring these as HHM(n), the bulk S'-matrix of the 
HHM(n) can be obtained from that of the RSOS(2n+l) by folding the heights as 

a, 2n + 1 — a — > n — a ;1 < a < n . (4-10) 
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The above one is well-defined without ambiguity, due to the symmetry of the S'-matrix. 
The integrability conditions such as the YBE and BYBE are transformed accordingly 
maintaing the structure. These mean that we can write the reflection amplitude of 
HHM(n) as 



R a bc (8) = 11(6) ( 



PbPc 



-6 /2m 



\PaPa, 

[6 b ^ c X£ c (u) + 5 bc {5 b>a+l U a {u) + 5 b>a V{u) + 6 bia ^D a (u)}] , (4.11) 



where p a denotes the q- number, 



sinaA _ , tt 

Pa = —. — r - a = n — a 



sin A 2n + 1 

Then the BYBE's solution for HHM(n) are 



U a (u) = (-l) a +Msin(2u-aA) + 



sin 2u 

D f sin(2^-aA) _ a 
sinaA \ sin2u 

D a {u) = (-lf +1 ,4sin(2u + aA) + 



+ D rsin(2u + aA)_ ( _ i) 



sin 2u 

iA) _ 

sinaA t sin2u 

V(u) = (-l) n+1 Asm(2u + n\) + —^— 

sin 2u 

| D f sin(2u + nX) _ 
sinnA 1 sin2rt 



x a-i,a+i x a+i,a-i = A 2 {sin 2 A - sin 2 aXj - 2AB {cos A - (-if +1 cos aXj 

„„ 9 f 1 11 „ 9 f cos A . cosaA 

+2D 2 — 2~r - — 2"~r \ + 2D 2 \ — 2~r - (-l^-^— 
{ sin A sin aX ) I sin A sin aX 

which can be read from eqns.( ^19| ), (|2.20 ) with odd a. 



5 Conclusion 

In this paper we derived the most general nondiagonal reflection matrices of the RSOS/SOS 
models and the hard hexagon model using the boundary Yang-Baxter equations. We 
find new one-parameter family of the reflection matrices for the RSOS model which gen- 
eralizes the previous result in H where there is no free parameter. This free parameter 
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can be used to control the flow between the fixed and free boundary conditions. Since 
the bulk RSOS theory describes the perturbed conformal theories by the least relevent 
operator, the boundary conditions can show how the conformal boundary conditions 
can change under renormalization group flows. 

It is still open problem to show the vertex-face correspondence under the presence 
of the boundary between the sine-Gordon and the SOS theories while our three classes 
of the SOS reflection matrices may be useful for this purpose. 

For the hard-hexagon model which can be mapped to RSOS(5) model by folding four 
RSOS heights into two, the solutions can be obtained similarly with a main difference 
in the boundary unitarity conditions. Due to this, the reflection matrices can have two 
free parameters. We show that these extra terms can be identified with the 'decorated' 
solutions. This means the general solution space of the BYBE is spanned by each 
fundamental solutions and decorated ones. 

Considering that the HHM is related to the perturbed conformal theory by the most 
relevent operator, it will be interesting to consider how the two different perturbations 
can make difference in the boundary interactions. 
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